INSTTTJ 
AFOSR-TN-59-944 


NEW  YORK  OiM 


25  Waverly  FLcc,  New  York  3,  N.  t 

^^EET  PR/J. 
/>       ^^   Jv;^     NEW    YORK    UNIVERSITY 

a:    I    iJ   Pf     ^     Institute  of  Mathematical  Sciences 


m 


^       ^        .^      Division  of  Electromagnetic  Research 


RESEARCH    REPORT    No.    BR-31 


On  Some  Fredholm  Integral  Equations 
Arising  in  Diffraction  Theory 


C.  H.  YANG 


Mathematics     Division 

Air    Force    Office    of    Scientific    Research 

Contract    No.    AF    49(638)229 

Project    No.    47500 

AUGUST,    1959 


^ 


HEW  YORK  UNTVERSITT 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 
Jlesearch  Report  No.  BR-51 


On  Some  Fredholm  Integral  Equations 
Arising  in  Diffraction  Theory 
C.E.  Yang 


C.H.  Yang 

aey~Bo^owitz  / 

Acting  Project  Direct*: 

August,    1959 


Qualified  requestors  may  ohtain  copies  of  this  report  from  the 
Armed  Services   Technical  Information  Agency,    Arlington  Hall 
Station,    Arlin^on  12,   Virginia.      Department  of  Defense   con- 
tractors must  be  established  for  ASTIA  services,    or  have  their 
'need  to  know'    certified  by  the   cognizant  mllitaiy  agency  of 
their  project. 


1. 


Abstract 


Kernels   of  integral  equations   arising  in  diffraction 
theory  are   investigated.      These  kernels  depend  on  a  parameter 
a.      Fredh.olm's  theory  shovs  that  the  resolvents   of  the  kernels 
are  meromorphlc  functions  of  a  and,   hence,    can  be  expanded  into 
power  series   of  a  vhich  are   convergent  for  certain  regions. 
Regions  of  convergence   for  a,    in  which  the  perturbation  method 
would  work,   are  examined. 
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2. 


In  order  to  show  thflt^  in  principle,  it  might  be  that  oiir  Fredholm  equa- 
tions have  a  unique  solution  for  all  a,  we  give  an  example  of  an  integral 
equation  of  the  type  studied  here  for  which  this  is  actually  the  case. 
(See  k.   Remarks  and  Examples,  P-  H) 


2.   Integral  Equations 

Let  us  consider  the  following  integral  equations  which  were  derived 
by  Magnus  L-^J  and  Jones'-  ■' ,  respectively. 

1 
(1)       S(x)  +  J    G(x,y;a)  S(y)  dy  =  T(x) 


where 


nr  \  -^      r     h  2      I    ^locwlx-yl        -iccw(x+y)    I     dw     . 

G(x,yja)    =       —     //1-w       <e        i-'i-e        ^      "    \     —     \ 


(2)  S(x)   +     r    K(x,yja)   S(y)   dy  =  F(x)    ,   S(-x)   =  -  S(x) 


-1 


where 


S(x)  is  vmknowQ  and  T(x),  F(x)  are  given  functions. 

In  these  integral  equations,  we  find  that  the  kernels  contain  a  para- 
meter a  in  a  non-linear  manner.  It  should  also  be  noted  that  both  kernels 
are  symmetric  in  x  and  y  but  that  they  are  neither  hermitian  symmetric  nor 
normal,  except  for  special  values  of  a. 

It  is  known  from  the  Jredholm  theory  that  (l)  and  (2)  have  the  unique 


3. 


2  2 

solution  S(x)  €  C(or  L  )  if  T(x),F(x)  C  C(or  L  )  and  if  the  homogeneous 

equations  have  no  solution  S(x)  €  C(or  L  )  except  the  trivial  solution, 

2 
where  C  and  L  are  classes  of  functions  which  are  continuous  and  square 

summable,  respectively. 

In  other  words^  if  a  is  not  an  eigenvalue  of  one  of  the  integral 

equations,  there  exists  a  unique  solution  for  that  equation.   If  this  is 

the  case,  then  the  \mlque  solutions  are,  respectively: 


J- 
(5)       S(x)  =  T(x)  -  J   R^(x,y;a)  T(y) 


iy 


W 


J. 
S(x)  =  F(x)  -   /  I^(x,y;a)  F(y)  dy 


where  R(x,y;a)  is  the  resolvent  for  the  respective  kernels. 
The  following  results  are  obtained. 

[51 
The  case  of  Magnus  '-^-' ; 


There  is  no  eigenvalue 
Im  a  >  0 
Re  a  =  0 


for 


V 


kl  <  % 


a  can  be  taken  as  1.24,  but  the 
actual  least  upper  bound  for  la| 
is  unknown. 


h. 


The  case  of  Jones 


[2] 


There   is  no  eigenvalue 
Im  a  =  0 
for       s        Re  a  =  0  and  b  >  0 


|a|      <  a-^  , 


where  cu,  >  1.19^  hut  the  least 
upper  bound  for  |a|  is  unknown. 


a-plone 


5.   Proof  of  the  Results  Obtained 

[5I  1 

The  case  of  Magnus  '-^-'  :   Let  S  =  u  +  iv,  -  G  =  K  +  iL,  a  =  a  +  ib 

a 

1 
where  u,v,K,L^    a  and  b   axe  all  real.     Let  us  also  define  Ef  =      /   H(x^y)   f(y)   dy, 

1  o 

(f^g)  =  /  f"(x)  g(x)   i±x.,   where  the  bar  indicates  the  complex  conjugate.   Then 

o 
we  have,  from  (l),  by  putting  T  =  0, 

(S,S)  +  (a  +  ib)  f (K  +  iL)S,SJ  =  0  . 

Separating  the  real  and  imaginary  parts _,  we  obtain 


(5) 


(S,S)  +  a(KS,S)  -  b(LS,S)  =  0 
b(KS,S)  +  a(LS,S)  =  0  . 


Solving  for  (KS;,S)  and  (LS,S),  we  obtain  (for  |a|  +  O) 


(6) 


(KS,S)  =  -^  (S,S) 


(LS,S)  =  -^  (S,S) 


From  (l),   we   also  obtain 


L  =  L(x,y;a,b)    =  i    |  A^  je-^^^^-^y^cos   av(x+y)-e-^^l^-y ■  cos   av(x-y)  I  | 
o  ^  i 

K  =  K(x,y;a,b)   ='^   J,^I7  ie-^^^^^y^sin  aw(x+y). 
^  o  l 


-bw  x-y      .  I         1  V  d.w 

■e        '      '''sm  awl  x-y  I  ?  — 


If  we  define 


(8)  < 


T^(x, p)    = 


TgCx,  p)    = 


e      ^  cos  awx  ,    if  x  >  p 


0  ,    if  x  <  p 

e^"^P  sin  awx  ,    if  x  >  p 


0 


,    if  X  <   p 


We  now  obtain 


(9)     L(x,yja,b)    =  -    ryi-w^  e"^''^''"'^^  Id  J  (t^(x,  p)T^(j,  p)^^{x,  p)l!^{y,  p))  dp 


+  —  sin  awx  sin  awy>  dw 


and 


where 


6. 


(10)  (LS,S)    =   (Lu,u)   +  (Lv,v) 


1  r      1 

o  L     o 


^+  f/)  \  dv 


(11)  < 


X 

r      0^   =  0^(x, p,w)    =      /     e"  "^  T^(x, p)   u(x)   dx 

o 

1 

f.    =  T.(x,p,v)    =      r    e"^"™  T^(x,p)  v(x)   da 

o 

1 

F     =  F  (x,v)    =      /   e  sin  avx  u(x)   dx 


(i   =  1,2) 


L      Fg   =  F2(x,v)    =      /   e"  ^^  sin  avx  v(x)   dx      • 


Hence   (LS,S)  <  0  If  b  >  0;  but  from  (6)^   b  >  0  implies   (LS^S)   >  0  and 

It  follows  that   (LS,S)    =  0  if  b  >  0.      Therefore,    (S,S)    =  0  if  b  >  0,   by  (6) . 

If  b   =  0,   we  obtain,    from  (5),    (LS,S)    =  0  (a  #  0,    otherwise  trivial)  which 

lirrplies  F  =  Fp=  0  by  (lO).      It  follows  that  u  =  v  =  0. 

If  a  =  0,   we  obtain   (S,S)    =  b(LS,S)   by  (5)    and,   therefore, 

1  1 

(S,S)   +  b^      \K^     r^^A  *2^)   dp  dw  =  G 


by  (10)    and  (ll).      Since  both  terms   on  the  left-hand  side   are  non-negative, 
we  obtain  (S,S)    =  0. 


7. 


1  1 


It  is  known  that  if  /  /  |G(x,y;a)  j^dxdy  <  1,  then  the  Reumann  series 
o  o 
of  (l)  is  convergent;  that  Is^  for  a  satisfying  the  above  condition,  R  is 

regular  and  a  is  not  an  eigenvalue  of  (l).   Since 

2      ^  (x+y) 

G(x,y;a)       "^   J    JI^    J  e^^P  dp  , 

o  ix-yi 


|G(x,y;a)l      =  G(x,yja)    G(x;,y}a) 


i^       1     (x+y)  1     (x+y) 

i^    J  JZ^   J  e^«^P  dp    J   JI^   J  e-i^^^  dX 

o  |x-y|  o  lx-y| 


1,     /     1    , (x+y)  \i 

<14.^vC7dv|  e-^-Pdp) 


-  (4b)- 


therefore^ 


11  1^         1  1 


//     lG(x.y;a)l2d.dy<     -^//    (e  1^  I  C^^)-   e  1^  l-H  f^  .^ 


(MH) 


<     -^ 


(M^i) 


J    e^l^I    +4(l-2|b|)e^i^I    -  l^|b|    -  5    I 


=     f(a,|b|)    . 


8. 


By  putting  f(a^ lb])  <  1,  we  find  that 

vhen  b  =  0  ,  a  =  2.23 
vhen  a  =  0  ,  b  =  1.24   . 

[2] 
The  case  of  Jones  *-  -^  :  Let  S  =  u  +  Iv,  K  =  A  +  IB,  a  =  a  +  ib^  where  u,v,A^B^ 

a  and  b  are  all  real  as  before.   Then  from  (2),    F  =  0^  (S,S)+((A+1B)S, S)  =  0; 

separating  the  real  aad  imaginary  parts ^  we  obtain 


(12) 


(S,S)  +  (AS,S)  =  0 
(BS,S)  =  0 


where  A  and  B  are  as   follows; 


X 

A  =  A(x,yja^b)   =  —  cosh  a(x-y)  ^ — ^  =  p—  cosh  a(x-y)/    cob  b(x-y)tdt 


(15)      { 


-„       T,/              -u  \        -1     sinh  a(x-y)  ,  /        \ 

B  =  B(x,y;a,b)    =  —    ^_^    ^'    cos  b(x-y) 


(13)   is  derived  from 


K  =  A  +  iB  = 


_^     ga(x-y)   ^ib(x-y)_  ^-a(x-y)   ^-ib(x-y) 


2rt 


x-y 


If  b  =  Im  a  =  0,  theo  from  (15),  we  have  A  =  0.  Hence,  (S,S)  =  0  by  (12). 
There  is  no  solution  except  the  trivial  solution  for  Im  a  =  0.   If  a  =  Re  a  =  0^ 
then 


A  =  |-  r  cos  b(x-y)tdt  , 
-1 


by  (15).      It   follovs   that 


(AS,S)    =   (Au,u)    +   (Av,v) 


1 

=  h  J  (\'*  '2'*  '1* ''/'  ^^ 

-1 
vhere 

1  1 

0^  =   /  cos  bxt  u(x)  dx  ^   Op  =   /  sin  bxt  u(x)  ix 

-1  -1 

1  1 

1     =      /  cos  bxt  v(x)  dx  ,   ^  =   /  sin  bxt  v(x)  dx  . 
-1  -1 

Hence,  (AS,S)  >  0  If  b  >  0,  in  which  case  ve  obtain  (S,S)  =  0  from  (12), 
From 


I  ,    ,|2  _   1   slnh  a(x-y)    slnh  a(x-y) 
mx,y)|   -  ^  —-  ^ 


-y  x-y 


^rt^(x-y) 


X-V^  N  '^  ''^ 


^ 2    f cosh  2a(x-y)  -  cos  2b(x-y)  j 

'x-y)     V  / 


2jt"(x-y) 


J:^   y  ^2a)2°+  (-l)^^l(2b)2"  V  li^=^^ 
2-   4  I  J     (2n)' 


2n-2 


10. 


ve   obtain: 
1     1 


■1    '-1  '^'^       n=l      L  J 


1     1 


X      /      /      (x-y)^""^  dxdy 


-1     -1 


2n+l 


= ;?   I  j'"'"*  (-ir^(-)^l  pirr  lln 


^ 


<     f(a,b) 


vhere 


f(a,b)    =  ^ 


Yp   (cosh  ha  +  cos  ifb)    +  —  (a  -b    )    +  Vc 


X 


(,„,,, ,.x-i||)!-ijt^)] 


Putting  f(a,b)  <  1^  we  obtain  a  <  1.19  when  b  =  0  and  b  <  1.5  when  a  =  0. 
If  we  take 


f(a,b)  =  4  (  ^  =°^^  ^a  -  I  +  ^^^)   , 


-(  - 

2  ^  2 


then  we  have  ^  <  ^  when  a  =  0. 


11. 


h.        Remarks  and  Exanrples 

We  may  summarize  the  methods  used  here  as  f ollovs :  Let  us  consider 
the  following  Integral  equation 

(Ik)  S(x)  +  j^    G(x,y;a)  S(y)  dy  =  T(x) 

o 

■with  the  kernel  G(x,y;a)  =  2a  e""^^"^^^  . 

We  can  easily  find  the  Fredholm  determinant.  Its  First  Minor  and  the 
resolvent  of  (l4)  as  follows: 

D(a)  =  e^"- 

D(x,y;a)  =  2a  e^^^^^^^ 

vfv  -^-rr^    -   P(^^y^°')  -  p„  a(x+y-2)  , 

In  this  case,  the  Integral  equation  does  not  have  any  eigenvalue  for 

any  finite  a  since  the  Fredholm  determinant  D(a)  is  an  entire  function 

which  has  no  zeros  for  any  finite  a,    and  we  have  the  unique  solution: 

1 


S(x)  =  T(X)  -   r2a  e^'^^'-y-^)  T(y)  dy  . 


In  the  example  above,  it  is  easy  to  compute  the  Fredholm  determinant: 
however,  we  are  dealing  here  with  a  rather  special  case.   In  mxjst  cases 
(see  Magnus  and  Jones),  it  is  very  difficult  to  evaluate  the  Fredholm 
determinant.   Since  we  know  that  zeros  of  the  Fredholm  determinant  are 
eigenvalues  of  the  integral  equation  from  which  it  arises,  one  method  for 
locating  eigenvalues  is  to  Investigate  zeros  of  the  Fredholm  determinant. 


12. 


If  an  evaluation  of  Fredholm's  determinant  Is  difficult,  an  indirect 
method  proving  non-existence  of  eigenvalues  for  certain  regions  is  useful. 
If  the  real  or  imaginary  part  of  a  kernel  considered  is  positive  (or  nega- 
tive) definite,  the  indirect  method  would  -work  out  in  general  (see  Magnus' 
case  for  Im  a  >  0,  Re  a  =  0;  Jones'  case  for  Im  a  =  0  and  Re  a  =  0,  b  >  O) . 

Another  indirect  method  is  to  obtain  a  majorant  f(a,b) 

of  /  /  |K(x,y;a,b)|  dxdy  (see  Section  5)  an<3.  to  find  the  region 


R  =<   (a_,b)  |f(a,b)  <  1  >-  of  non-existence  of  eigenvalues.   The  region 
obtained  by  this  method  is,  hovever,  limited  to  the  neighborhood  of  the 

origin  and  varies  according  to  the  majorant  function  used. 

some 

It  should  also  be  noted  that  for^ kernels   of  the  type   considered  (symmetric 

but  non-linear  with  respect  to  the  parameter),    there  exist  no  eigenvalues  even 

if  ve  have   a  real  symmetric  kernel.      This   cannot  happen  if  K(x,y;a)=aX(x,y), ' 

where  K(x, y)   is   independent   of  a  since  we  know  that  there  then  exists  at  least 

one  real  eigenvalue  a      (e.g.    See    [6]). 

We  shall  now  consider  another  example  arising  from  the  work  of  Noble'-  -■  : 
1 
(15)  W(x)    +     fj  W(y)   K^(x,y;k)   dy  =  F(x)      (O  <  x  <  l) 

I 

k 
-i   j     H^^''(xt)j  (yt)  J^^tP  dt, 
o 
K^(x,y;k)    =  ^ 


where 


X  >  y 


-i   j     H^^^(yt)j  (xt)  /k^-  t^  dt. 


y  >  X 


L      o 


15. 


the  integral  equation  (I5)  is  one  of  Fredholm's  type^  and  symmetric  In 
X  and  y  with  the  weight  function  y. 

Let  W  =  u  +  iv^  K  =  S  +  iT  (for  k  real),  then  we  obtain  from  (I5), 
by  putting  F  =  0: 

r  u(x)  +  I  y  js^(x,y;k)u(y)  -  T^(5t,y;k)v(y)  j  dy  =  0 
o 


(16) 


v(x)  +   /  y  js^(x,y;k)v(y)  +  T^(x,y;k)u(y)  1  dy  =  0 
o 


where 


k 


T  (x,y;k)  =  -  /  /k^  -  t^  J  (xt)  J  (^)  dt 

M-  ,J  M-       M- 

O 


since 

H^"^Hxt)  =  J  (xt)  +  iW  (xt). 

After  eliminating  S  from  (16),  we  obtain 

(IT)      \     \     T^(x,y;k)  iyia(y)  xu(x)  +  yv(y)  xv(x)  >  dy  dx  =  0  . 

00  L  J 

Changing  the  order  of  integration^  we  obtain 

k    r      1  2  1  2] 

r  /k^  - 1^  <  ( rxu(x)j  (xt)dx^  +  ( rxv(x)j^(xt)dx^  idt  =  o 


Hence 


(18)    / 


/  xu(x)j  (xt)dx  =  0 
o 

/  xv(x)j  (xt)dx  =  0  , 


"^      o 


14. 


and  it  follows  that 


(19)  <       o 


X     ^         u(x)(ix  =  0 


X     ^         v(x)dx  =  0 


(n  =  0,1.2,...) 


We  obtain  u  =  v  =  0,   by  (l9),   vhenever  |a  is   real  and  |j.  >  -   2. 

When  \i  =  —  ,    (15)    can  be   reduced  to  Magnus'   equation  (l)  by  setting 
S(x)    =yxW(x).      In  this   case. 


Ki(x,y;k)    =  ^     — 

2  It        I — 

yxy 


r  ^^^"^^  /e^l^-y|t .  e^^^+y)*"^ 


dt 


1   A  t2   r 


dt 


=     G(x,y;k)/  y/^ 


15. 


Appendix  I 


Let  us   consider  the   follovlng  integral  equation: 
b 
(20)  S(x)    +  \  j        G(x,y;a)S(y)dy  =  T(x)  - 


By  Predholm's  theory,   we   obtain  the  determinant  D  (X),    the  First  Minor 
D   (x,y;X)    and  the  resolvent  R  (x^yjX)    as   follows: 


b       b 


n  =1  a        a 


a        a 
~n~ 


G(t^,t^;a),...,G(t^,t^ja) 


G(t   :,t,  ja),  .  ..,G(t   ,t   ;a) 
^   n     1  ^   n'    n 


dt^  ...dt 
1  n 


(22)      D^(x,y;X)    =  G(x,y;a) 


I  ^  /•• 


b       b 


n=l 


a        a 
n 


G(x,yja),G(x,t^;a),  . . .  ,G(x,t^ja) 

}  •  •  • }  •  •  • } 

} }  •  '  • } 

G(  t^y; a) ,  G( t^,  t^;a) ,  .  . . ,  G( t^^t^J c^) 


dt,  ...dt 
1  n 


(25)     Rjx.yjX)    =  Djx,y;X)/DjX) 


If  X  =  -  1,  then  (20)  is  of  the  type  considered  here.   Therefore,  ve 
define  D(a),  D(x,y;a)  and  R(x,yja)  by  the  identities: 


(2M      D(a)  H  D  (-1) 
a 

(25)  D(x,y;a)  ^  Djx,y;-l) 

(26)  R(x,y;a)  h  R^(x,yj-l) 


16. 


If  G(x,yja)  is  an  entire  function  of  a,  then  D(a),  D(x^y;a)  are  entire 
functions  of  a  and  R(x,y;a)  is  a  meromorphic  f-unction  of  a. 

Since  D(a)  is  an  entire  function,  Picard's  theorem  shows  that  D(a) 
has  infinitely  many  zeros  unless  zero  is  the  exceptional  value.   Therefore, 
except  in  this  case,  R(x,y;a)  has  infinitely  many  poles  which  are  eigen- 
values of  the  equation  (20)  with  X  =  -1. 

Since  D(a),  D(x,y;a)  are  entire  functions,  D(x,y;a)  can  he  expanded 

into  a  power  series  in  a  which  converges  uniformly  and  absolutely  for  all 

finite  a;   l/D(a)  can  be  expanded  into  a  power  series  in  a  which  converges 

uniformly  and  absolutely  for  jal  <  Min  |a. [  where  a.  are  the  zeros  of  D(a) 

i 
(i.e.,  eigenvalues). 

Hence,  R(x,yja)  can  be  expanded  into  a  power  series  in  a  which  is 

uniformly  and  absolutely  convergent  for  |a|  <  Min  ja. |  . 


See,  e.g.,  G.  Vallron  -  Theory  of  Integral  Functions; 
New  York,  Chelsea  Publ.  Co. 


17. 


Appendix  II 

For  the  case  of  diffraction  by  a  circular  aperture,  the  eigenvalues 
and  eigenfunctions  have  been  calculated  approximately  by  means  of  the 
geometrical  theory  of  diffraction.   This  method  shows  that  the  eigenvalues 
are  asymptotically  equal  to  the  roots  of  the  follovring  equations  (See  [5], 
pp.  1+55-458) 


(27)      (2.a)'  e-^(2a  +  n/h)^  ^   ^  ^_ 


Equation   (27)   has   infinitely  many  roots  which  are 
a  =  -    .l^i   ,    -    Ali 

i(   —  +    (?    ) 

a=-r-e  (O^-asn-*") 

n        4  ^      n       2 


where 


(5-     +  2n«   =  7     sin    Q- 


log  jt   7        =7     cos    Cj    . 
'n  n  n 


i 

b 

a-plar 

• 

• 

• 

• 
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